Thermal transport properties of a charge density wave 
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Effects of collective modes on thermoelectric properties of a charge density system are studied. 
We derive the temperature dependence of thermoelectric power and thermal conductivity by apply- 
ing the linear response theory to Frohlich Hamiltonian. Energy dissipation has been attributed to 
nonlinear interaction between phase mode and amplitude mode, ignoring disorder effects. We have 
found that the temperature dependence of the correlation function of electrical and heat currents 
is the same as that of electrical conductivity. This implies that thermoelectric power is inversely 
proportional to temperature. We have also found that temperature dependence of all the correla- 
tion functions are essentially determined by the common mechanism - nonlinear amplitude-phase 
interaction. Thermal conductivity has nearly constant value at the temperature above amplitude 
mode gap, and has exponentially low value at the temperature sufficiently below it. 

PACS numbers: 71.45.Lr, 65.90.-K, 05.60.Gg 



I. INTRODUCTION 

For many years, thermoelectric properties of various 
conductors have attracted both theoretical and experi- 
mental interest. These properties have also been inves- 
tigated from practical standpoints. In this study, we 
are interested in thermoelectric transport properties of 
a charge density wave (CDW) system especially in con- 
nection with possible effects of collective modes. 

As for electrical conductivity, a lot of studies on CDW 
system have been done until now. One of the dominant 
interest of such studies are focused on collective-mode 
contribution to conductivity. CDW has two collective 
modes- amplitude and phase modes. As is well known, 
phase mode can carry permanent currents in principle 
0,0- In practice, however, impurity and/or commen- 
surability pinning prevent this phenomenon. Some ef- 
fects of phase mode have been found in experiments 
of nonlinear conductivity and frequency dependence of 
ac-conductivity. Moreover, certain features of electrical 
dc-conductivity in CDW systems could be explained as 
many-body effects arising from nonlinear interaction be- 
tween phase and amplitude modes 

In contrast to electrical conductivity, theoretical study 
of thermal transport of this system has been rather rare. 
Experimentally, evidence of collective mode current of 
thermoelectric power have been found on nonlinear re- 
sponse to electric field 0, 0|. Some experimental data 
have also suggested fluctuation effect in thermal conduc- 
tivity near transition temperature 0,0,13- At present, 
it is not easy to observe collective-mode contribution in 
thermoelectric response, especially at low temperatures. 
This is because of pinning effect which is seen in electrical 
conductivity as well. However, in principle, these effect 
should exist in this region, and indirect evidences are ex- 
pected to be actually observed. Theoretical study of this 



region is expected for better understanding of thermo- 
electric transports caused by collective modes. For these 
reasons, we investigate the collective mode response to 
thermoelectric power and thermal conductivity of this 
system by applying linear response theory to a micro- 
scopic model. 



II. BASIC FORMALISM 



We start from the Frohlich Hamiltonian 



^ = E« 



CfeCfcCfc + w q b\b q 



■bi q )cl +q c k . 



(1) 



qk 



Here, electron-phonon interaction coefficient 7 are as- 
sumed to be constant equal to the values at \q\ = Q = 
2kp, since only the small portion of the phonon state 
are relevant in considering the CDW collective modes. 
We ingore electron spin indices for simplicity. Adopting 
Nambu representation^ @> we can rewrite the above 
Hamiltonian as follows, 

H = J2cHm k r 3 + ^-r )c(k) + nJ2 E 



Q 3=1,2 



+7feE E h(q)+b](-q)]c\k + q )T jC (k),(2) 

V iN k ,q .3=1,2 

where 7i,T2, T3 are Pauli matrices, tq is the unit matrix, 
6f is Fermi energy and 
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Summations are done in the range of |fc| < Q/2 and 
| q | < Q with respect to electron and phonon states. 

The term is comes from the quadratic dispersion. 
Since transport phenomena are usually related to states 
near the Fermi surface, linear dispersion relation is of- 
ten adopted, espesially in one dimensional models. How- 
ever, as far as thermoelectric transports are concerned, 
electron- hole asymmetry turns out to be important. It is 
because particles and holes respond to electric field in an 
opposite way: total energy flow vanishes if we have the 
exact electron-hole symmetry . 

Electrical conductivity a is given as follows 



T 



(4) 



where T is temperature and L 11 is transport coefficient 
calculated from current-current correlation function us- 
ing Kubo- formula |ld| | . Similarly, thermoelectric power S 
and thermal conductivity k are given by 



S 



J>2 1 



eT L 11 ' 



(5) 
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where L 12 is calculated from correlation function between 
current and heat current operator, and L 22 correlation 
function between heat current operators, |l2|.and 
[l3| . for example. 

Electric current operator is given by 



J = ev F S2 c f (k)(T 3 + ^To)c(fc). 



(7) 



On the other hand, the heat current operator is composed 
of two parts: electron part and phonon part. Histori- 
cally, there are several forms that represent phonon part 
[l4| . Such difference originates from various treatments 
of electron-phonon interaction. The following discussion 
is hardly affected by such differences. In this study, we 
give heat current operator of electron part as 



J Q - 

v F i(e n + uj l /2)c){k, e n )(r 3 + -^-T )c(k, e Tl 

, — 



•wj). (8) 



where, e n and ui n are Matsubara frequencies for fermions 
and bosons. Note that, we need not to define heat current 
operator as Jq — fiJ since we have already chosen the 
Fermi energy as the origen of energy. Phonon part is 
also given by 

Jq 1 = v ph ^2 *( £ ™ + ^/ 2 ) 

<7>0,n 

x [b\{e n )b q {e n + u{) - &t 9 (e n )6_ g (e„ + 0Ji)], (9) 
which is rewritten in 

Jq 1 = I'ph *( £ " + W // 2 )^I(9, e n)6j(g, £n (10) 



where, both e„, and toi are boson Matsubara frequencies. 
Using above operators, transport coefficient, i IJ is given 
after usual analyitic continuation to real frequencies, by 



L lJ = lim Im 



L y '(w + iOJ 



(11) 



where L lJ (io>„) is corresponding correlation funciton. 

Here, We restrict ourselves to low temperature region 
T « T p in the following discussion. Peierls gap A de- 
fined by 



A = 



1_ 

N 



(12) 



is well-developed in this temperature range. Under such 
conditions, b\(q) and 62(9) can be identified with the an- 
nihilation operators of the amplitude and phase modes 
of CDW. Electron Green's function in the mean field ap- 
proximation is 



G(p,ie n ) = 



1 



tie - -^) 2 - f 2 
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Collective modes are described by the following Green's 
functions 



D a (q,iu n ) 
D^q, iu) n ) 
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-<T T a g (T)a- q {0)>, (14) 
T <T T <p q (r)(f>- q (0)>, (15) 



here, ot q (r) and <pq are displacement operators of ampli- 
tude and phase modes given by 



u q = bi(q) + b[(-q), 
4> q =b 2 {q)+b\{-q). 



(16) 
(17) 



In the low energy region, above Green's function can be 
written in 



D a (q,iuj n ) 



2Q 



2fi 



D^{q,iu} n ) 



(18) 
(19) 



where Q is co±2k F i A is dimensionless coupling constant 
7 a a is lattice constant, and /1 = i^r is the CDW 



n-Kv F ' iD WJ """' 1 ' i f - AO- 

mass paramter. Ordinarily it has quite large values typ 

ically on the order of 10 2 . 



III. THERMOELECTRIC POWER 

Now let us calculate thermoelectric power. In this 
study, we have a interest to collective mode response, so 
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FIG. 1: Feynman diagrams which give phason self-energy. 
Each solid and dashed lines show contribution of phase mode 
and amplitude mode. Electron vertex is shown as triangle and 
square points whose contributions result in constant values. 



we don't consider impurity effects. We also neglect effect 
of quasi-particles. Needless to say, these have much ef- 
fect in real systems, but at least, these contributions are 
expected to be vanished at low temperatures T « Xp. 
Damping is derived from nonlinear interaction between 
phase mode and amplitude mode, which had been devel- 
oped by one of the author @, || , yielding the following 
expression for the phason life time from Fig. 1 



^#Tsh 2 (^), 



There is a factor in above two equations. This means 
that L 12 and L 22 can be non-vanishing only when the 
electron- hole symmetry is violated. It is derived from 
asymmetry between electric and heat current operators. 
Note that, the factor ^ is not necessarily small in CDW 
system : a typical value is on the order of 0.1. 

Using Eq. (4), (5) and 1)21(1 . see that the thermoelectric 
power is given by 



2e X'e F T 



(23) 



(20) 



In the experiment on TaS3, field dependence of thermo- 
electric power corresponded to differential resistivity Q|. 
This suggest that L 12 doesn't have dependence on collec- 
tive modes. However, in the system of the blue bronze, 
nonlinear field dependence of L 12 has observed at sev- 
eral temperatures 0. In this system, temperature de- 
pendence of thermoelectric power which is shown in Eq. 
(23) is expected to be observed under the high electric 
field: above the threshold. We also note that Eq. 1(23(1 
diverges at zero temperature, which may be an artifact 
of our assumption: we have ignored impurity effects and 
quasi particle contributions. 



where, ui a = y 4^fi. Eq. ((20(1 shows that r is inversely 
proportional to T at temperatures above amplitude gap 
uj a , and diverge exponentially in the zero temperature 
limit. 




FIG. 2: Feynman diagrams which give transports coefficients. 
Double lines represent electrons. One-loop-type diagram (a) 
contributes to electric conductivity and thermoelectric power, 
while Aslamazov-Larkin diagrams (b) ~ (e) contribute to elec- 
tron part of thermal conductivity. 

In the previous paper electrical conductivity a is 
calculated from the diagram shown in Fig. 2 (a), and is 

given by o — jf^T. This diagram also contributes to 
other transport coefficients L 12 and L 22 . We denote these 
contribution as L\ 2 and L 22 which are given by 
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IV. THERMAL CONDUCTIVITY 

Thermal conductivity k is composed of two parts. We 
denote the contribution of electron and that of 

phonon as « p h 

First, we investigate the electron part. We see that 
diagram shown in Fig. 2 (a) doesn't contribute to collec- 
tive mode response of thermal conductivity. Actually 
we can easily confirm the relation L\ 2 — (L\ 2 ) 2 / L 11 = 
from Eq. 1(21(1 and Eq. ((22(1 . This represents that phase 
mode transport which is damped by nonlinear interaction 
between phase mode and amplitude mode is completely 
canceled out by back flow, i.e. contribution of thermo- 
electric power and electric current, which is shown in Eq. 
(6). 

Thus, we investigate other processes of the thermal 
conductivity. As for higher order diagrams, we focus 
on Aslamazov-Larkin-type diagrams. Note that Maki- 
Thompson and density of state diagrams are both im- 
portant factor in considering the superconducting fluc- 
tuations. But these terms give no contribution in our 
case. We find that Fig. 2 (b) ~ (e) give finite values to 
thermal conductivity, and the other diagrams vanish. As 
for electrical conductivity and thermoelectric power, all 
diagrams in this order don't have finite value. 

Actual calculation of these diagrams is rather compli- 
cated. We calculate these terms by the following way. 
Summation of these terms are represented as 

/( w FQ , ,iw«)[£>< ? i(wF'5 , ,i^n)£> ce (w_F9,i^n + m) 



(24) 
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FIG. 3: An example of diagrams which are shown in FIG. 2. 
Suffix k and q are momentum, ie n , iw m and iui are Matsubara 
frequency. is analytically continuated to lo + iO. 

where N is number of lattice site, and — 1/T, function 
f(q, iio n ) is contribution of electron parts shown in Fig. 
3. As is shown in Fig. 3, we have to consider two summa- 
tions over Matsubara frequencies ie„ and iu) m . However, 
we only pick up poles of Green's functions of phonon with 
respect to iu> m , since we consider the temperature suffi- 
ciently below the CDW gap, and effect of quasi particle 
excitation is expected to be very small. In this assump- 
tion, transport coefficient derived from Eq. H24|) is given 
by 

"^T^mlV^ 

Function f(\/6A,uj a ) is still hard to calculate directly. 

However, we find /(0, w Q ) + /(0, -w„) ~ ~(|^) 2 ^(1 + 
o(^)). Thus we use this value and finally find to be 

Next, we calculate phonon part. We find Eq. I|1U|) 
represents decay process between phase and amplitude 
mode, and it's correlation function is composed of prod- 
uct of Green's function of phase mode and amplitude 
mode. We calculate one-loop contribution of thermal 
conductivity of phonon part which is given by 

k c \ is similar to k p ^, both of them are inversely pro- 
portional to phason life time r. Actually, we substitute 
f(vFq, iu) n /2) with (\ui n -l-iu;;) 2 , then K p h can be derived. 
We may consider that these heat currents are essentially 
same: they are carried by energy-interchange process be- 
tween phase and amplitude modes. Such interpretation 
is reasonable, because the heat transfer is process of en- 
tropy diffusion. If local amplitude mode excitation causes 



phase mode excitation and, thereafter phase mode bring 
about amplitude fluctuation in some other place, it is 
interpreted as heat current. 

Finally we estimate the relative importance of these 
contributions to thermal conductivity. k c \ and K p h has 
exponentially low value at temperatures sufficiently be- 
low amplitude gap while they have nearly constant value 
k c \ and K p h at temperatures above it. Ratio between 
these values and electrical conductivity for this temper- 
ature range is given by 

1 / 3A \ 2 

where Lo is Lorentz number -^-(fce/e) 2 , fee is Boltzmann 
constant. Since a factor (A/cf) 2 is quite small which has 
an order of 10 -2 , electonic contribution is small. On the 
other hand, phonon contribution is much larger, because 
*x(v P h/vF) 2 is expected to have an order of 0.1, where we 
estimate v p h/vp have an order of 10~ 2 . 



V. CONCLUSION 

In conclusion, we have studied thermoelectric proper- 
ties of a charge density wave system. We have particu- 
larly focused on collective mode response. We have found 
temperature dependence of correlation function between 
electrical current and heat current is same as electri- 
cal conductivity, which means thermoelectric power is 
inversely proportional to temperature. We have also 
found that thermal conductivity is calculated by energy- 
interchange process between phase mode and amplitude 
mode, and is inversely proportional to phason life time. 
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